ON FAST AND SLOW TIMES IN MODELS WITH DIFFUSION 
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The linear Kelvin- Voigt operator Ce is a typical example of wave operator Co perturbed by 
higher - order viscous terms as eu^xt. If Vs is a prefixed boundary - value problem for Ci;, when 
£ = Ci; turns into Cq and Vs into a problem Vo with the same initial - boundary conditions of 
Va- Boundary - layers are missing and the related control terms depending on the fast time are 
neglegible. In a sm,aU time - interval, the wave behavior is a realistic approximation of u^ when 
£ — > 0. On the contrary, when t is large, diffusion effects should prevail and the behavior of Ug 
for e — > and t -+ oo should be analyzed. For this, a suitable functional corrispondence between 
the Green functions Q^ and Qq of 'Pe and Vo is achieved and its asymptotic behavior is rigorously 
examined. For this, a suitable functional corrispondence between the Green functions 6e and Qq 
of Vn and Vo is derived and its asymptotic behavior is rigorously examined. As consequence, 
the interaction between diffusion effects and pure waves is evaluated by means of the slow tim,e 
et; the main results show that in time - intervals as (e, 1/e) pure waves are propagated nearly 
undisturbed, while damped oscillations predominate as from the instant t > 1/e. 

Green's Function, Partial differential equations, Viscoelastic models. Singular perturbations. 



1. Introduction 

Consider the class of dissipative phenomena described by the following model: 



£e-Ue = {edxxt + C^dxx - dtt)Ue = /, 



(1.1) 



where e, c are positive constants and the source term / may be linear or not. In the 
linear case, when a prefixed boundary-initial problem Ve is stated, the knowledge 
of the related Green function G^ allows to solve explicitly Vg . When the function 
/ is not linear, then G^ represents the explicit kernel of the integral equation to 
which the problem Vg can be reduced. 
Two typical examples in the non linear case are: 

i) For / = auf + bsinu, the equation (1.1) is the perturbed Sine-Gordon 
equation which models the Josephson effect in Superconductivity [1]. 

ii) For / = f{ux,UxtTUxx)^i^s equation (1.1) is the Navier-Stokes equation 
for a compressible gas with small viscosity [2] . 

Moreover, as for the artificial viscosity methods, (1.1) represents a model of 
wave equations perturbed by viscous terms with a small parameter e [3], [4]. In 
this framework, the behaviour of Me as e — > must be examined and the interaction 
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of pure waves with the diffusion effects caused by euxxt must be estimated. This 
interaction is meaningful in the evolution of many dissipative models (viscoelastic 
liquids or solids [5] - [9], real gases with viscosity [10], magnetohydrodynamic fluids 

[11]-) 

For £ = 0, the parabolic equation (1.1) turns into the wave equation 

A^o = {c^dxx - dtt)uQ = /, (1.2) 

and Ve changes into a problem Po foi' Wo, with the same initial-boundary conditions 
of Ve- So, boundary-layers are missing, but the approximation of Ue by mq is rough, 
because in large time-intervals the diffusion effects are dominant. As consequence, 
the asymptotic behaviour of u^ should depend on the slow time e t typical of the 
diffusion; on the contrary, boundary layer control terms characterized by the fast 
time t/e, should be negligible. 

Aim of the paper is to derive and analyze an appropriate functional relation 
between G^ and the wave Green function Go related to T'o (Th. 3.1). By means of 
this transformation, a rigorous asymptotic analysis of G^ is achieved (Th. 4.1), and 
Gg is approximated by solutions v{x, t) of the second-order diffusion-wave equations 

s s 

-Vxx = Vt+CVx, -Vxx = Vt-CVx (1.3) 

2 2 

which correspond to the heat equations Vyy = ve, where the time - variable is just 
the slow - time (9 = et) and the space- variables y = x±ct are related to traveling 
or backword waves. 

The physical meaning of the above analysis is clarified by the explicit solution 
related to the case / = (n.5). Except errors of order 0{e/t), this relationship is 

Ue{x,t)^ / e 2.t uq{x,t) dr, (1.4) 

V27ret J-oo 

and it clearly shows the interaction between diffusion and wave propagation. In 
the time - interval (£,-), when ei < 1 , pure waves propagate almost without 
perturbation; while as from the instant t > 1/e, damped oscillations predominate. 

Let us remark, generally, asymptotic theories are developed for models with first 
and second order operators as documented in [12] and [13]. For these operators, in 
fact, numerous maximum theorems for the rigorous estimate of the remainder are 
well known. 

As for the third - order model (1.1), the error of the approximation has been 
examinatcd by means of the functional relation (3.1). 



2. Statement of the problem 

If Ue{x,t) is a function defined in the strip 

Sl = {(x,t) ■.0<x<l, t>0}, 



let Ve the initial- boundary value problem related to the Eq. (1.1) with conditions 

f We(a;,0) == /o(x), atMe(x,0) = /i(x), xe[Q,l], 

{ 2.1 

[ Ue(0,i) = (^(t), uS,t)^i^{t), t>0, 

where /o, fi,ip,ip are arbitrary given data. 

Boundary conditions (2.1)2 represent only an example of the analysis we are 
going to apply; flux-boundary conditions or mixed boundary conditions could be 
considered too. Further, as f{x,t),fo{x),fi{x) are quite arbitrary, it is not restric- 
tive assuming ^3 = 0, V' = 0; otherwise, it suffices to put 

Ue^U^--1p —if] f=J + jV-\ ^if (2.2) 

and to modify /o(x),/i(x) consequently. 

Let z(s)dcnote the Laplace transform of the function z{t) and let 

F{x,s) = f{x,s)~sfo{x)~f,ix). (2.3) 

Further, if a{e) — s{es + (?Y^I'^ , let 

^(^'^)-2rT^f^%^ (2.4) 

lyes + &-)asenh{lu) 

and 

Ge(e, a;, s) - g{\x - CI, s) - g{\x + ^1, s). (2.5) 

Then, the Laplace tranform u^ of the solution of the problem V^ is given by 

(^ . „ . 

u,{x,s)^ / [F(e,s)+e/o(0]Ge(a;,C,s)de (2.6) 

Jo 

When £ = one has (t(0) = sjc and (2. 4), (2. 5) turn into 

. / X coah{{l - y)al(\ 

my, s) = (2.7) 

2, c s senn(is/c) 

Go(C, a:, s) = go{\x -^\,s)- go{\x + ^\,s). (2.8) 

Problem Vo is defined by (1.2), (2.1) and the Laplace transform iiQ of it's solution 
is: 

M^, s) = / F(e, s)Go(C, a:, s) de (2.9) 

By comparing (2. 4), (2. 5) with (2.7), (2.8), the following relationship between 
Gg and Go is pointed out: 

Ge(.) = -^ Go( ^^ )- (2.10) 



In order to estimate Gg (t) in terms of Go (t) and to achieve the rigorous asymp- 
totic behaviour of Ge{t) for e -^- , the corrispondence (2.10) must be analyzed and 
its inverse Laplace transform must be established. 

3. Functional dependence between Gs and Go 



If /q denotes the modified Bessel function of first kind and order zero, the inverse 
Laplace - transform of (2.10) is fully given by the following theorem: 

Theorem 3.1- The relationship between the Green functions G^{t) and Go(i) 
is: 



Geit)^ -=^ / e-""^ dv Go(u)/o( Vu{v - u) du (3.1) 

ey'iret Jo Jo e 

and their C- transforms satisfy (2.10) in the half- plane Re{s) > —(? je. 
Proof. For Reies + c^) > 0, one has: 



g-V- e2„2 g^ \/eS + c2« 



£( g-T^) = ^ (3.2) 



therefore the Laplace tranform of (3.1) is : 



G,{s) = J^ / e~^^^^^^'"dv I Go{u)Io{ — Vu{v-u)du (3.3) 

eves + c^ Jq 



^3 



2c2 



Go{u)du / g-f ("+«)v^3+^ /o( 7^) dv. 



sVss + c/ Jo Jo £ 

Considering that 



e 



-^Vss+c 



loi—V^) dv = - \ (3.4) 



e c Ves + c^ 



see ([14]), one has 



Gg(s) = / e V^^^ GoMdu (3.5) 



therefore the above formula coincides with (2.10). 



Consider now a first application of Theorem 3.1. By means of tfie transformation 
(3.1), it's possible to deduce explicit expressions for Ge from known formulae of Go- 
So, for instance, if Go is represented by its Fourier series 

o oo -I 

Go = — > — sin( — cnt) sin(— nx) sin(— nO, (3.6) 



CTT ^ — ' n 
n—1 



then (3.1) allows to make Fourier series of Ge explicit. It sufhccs to evaluate the 
transform (3.1) of the time-component Go,n oi Go'- 



Owing to 

1^1 sin 

if one puts 



Go,n{t) - sin(y c n t). (3.7) 



{ay)Io{b^y{2v-y))dy = 2sin{av) '^"^^^^L^^'^ (3.8) 

Va^ - 0^ 



a^TTcn/l, b = 2c^/e, k = 2cl/{T:e), LOn = v^l - (n/fc)^, (3.9) 

the integral transform (3.1) of Go,„ gives 

Ge.n[t) = e 21^ . (3.10) 

So, by (3.1), (3.6), (3.10), one obtains 

2 -^ -^„2 t . cnn 1 TT X . /TT <.N /0 11N 

Ge ~ — > e 2i2 giYii — :—uJnt) sm — nx sm(— nf . 3.11 

n—l 

and this formula represents the Fourier series of Ge previously established in [16] 
by Fourier method. 

4. Asymptotic behavior of Ge 



An additional consequence of (3.1) is the asymptotic analysis of Ge when e tends 
to zero. Referring to (3.1), let r — t/e and let 



riu,v,T) ^ ^^^ Goitu) /o(2cVV^)e-^'^[i+("+-)'/41. (4.1) 



^3 

/TT 

Then, the function Ge can be given the form 



Ge= I T{u,v,T)du dv, (4.2) 

where Q = {(u, u) e [0, oo] x [O.oo]}. 



So, Ge depends on e only through the fast time r. When e is vanishing and 
t > e, the parameter r goes to infinity and the asymptotic behaviour of integral 
(4.2) can be rigorously obtained by Laplace method [15] 

To establish precise estimates of the remainder terms, let 

x(t)=Xo t-1/3, a(r)-ao t-^'\ (4.3) 

where xo s-iid <7o are arbitrary real positive constants such that x ^^^d a are less 
then one for all t > 1. Further, let 

Oo - {(ti, i') e [1 - X, 1 + X] X [1 - a, 1 + a]} C Q. (4.4) 

The following Lemma holds: 

Lemma 4.1 - For all t > I, it results: 

\Ge-ll r{u,v,T)dudv\< ^ie-^^''''\ (4.5) 

where the constants A^, /x depend only on c,X0j (^o- 
Proof. For all real and positive z one has Io(z) < e^, so that by (4.1) one obtains 

|r| < -^ {cy^f e-"'^ '*("'") (4.6) 

with 

h{u,v)^l + {l/A){u + vf -2^/^. (4.7) 



Further, for all {u,v) e Q, one has 

hiu,v) = ^[("-1)' + {v- If + iV^-lfiV^+l)^ (4.8) 

+{V^-lf(V^+lf] > ^[{u - 1)2 + {v- 1)2]. 
Therefore, indicating by C a constant independent by r, one has: 

rdudv\<C f f e-'^[("-i)'+("-i)'l^ du dv. (4.9) 

Qo J JQ-Qo 

Hence, by standard computations, the estimate (4.5) is obtained. I 



The estimates of Lemma 4.1 allow us to define the following function 

/'-^ 2 / ,2 

e" " 2.? G'o(x,t) dT (4.10) 

-oo 



and to obtain the following result: 

Theorem 4.1 - For all {x,t) G 17, and t > e, one has: 

Ge{x,t)^H{x,t,e) (1 + pi) +P2, (4.11) 

where the remainder terms pi , p2 are decreasing functions of the fast time t — t/e 
such that 



\pi\ < fci r-i; p2 < k2 e-^''''\ (4.12) 

with the constants fci, ^2, A^ depending only on c,xo, <^o- 

Proof. By asymptotic formulae of Bessel functions, for real, positive and large 
z one has ([17]): 

/o(z)= ^L (1+ro) (|ro|<co z-i) (4.13) 

and so the integral of F on Qq (see (4.1), (4.2)) can be given the form 

Tdudv^— / Go{tu)du ^-^TTTT e-" ^ ''("^") di;, (4.14) 

27r Ji_^ Ji_„ (■uw)V4 

where h{u,v) is defined in (4.7) and the remainder term ri is such that |ri| < 
ci (r-y/ro)^^, according to (4.13). 

The function h(u,v) has a point of absolute minimum in (mq, ^o) = (1, 1), where 
it vanishes. Therefore, it results: h — Hq + h^ with 

h„ = ^[{u-lf + (v-lf]; \h,\<c4\u-l\ + \v-l\]\ (4.15) 

According to the Laplace method, [15], for large t , the dominant contribution 
is due to a neighbourhood of (uo, ^o) , so that it's enough to expand F near to this 
point. Further, if u and v are substituted by 

u— 1 = —= ui, v—l = —= vi, (4.16) 

the terms of third order like /i, are vanishing with order of t^^''^ ; besides, the 
limits of integration become (— oo, +cx}) as: xV^ ~ Xo t^ and Cy/r = ctq t^'^ . 
Then, for t — > oo, one has: 

Vdudv^—l e-'^''^Go{t + Veiu)du e'"^"^ dv + 0{-). (4.17) 

The above formula, together with the results of Lemma 4.1 (see(4.5),(4.10)), implies 
(4.1) and the estimate (4.12)2 for p2- Finally, by means of routine calculations, the 
order of the error pi can be rigorously specified according to (4.12)i. I 



5. Diffusion and waves 



To remark the physical meaning of the results, let consider the simple case: 
/ = 0, /o = 0, /i 7^ 0. By (2.6) and (2.9) one has: 

ue^- f fi{OGe{tx,t)d^; uo^- [ /i(OGo(e,x,t)dC (5.1) 

Jo Jo 

Consequently, except errors of order 0{e/t), Theorem 4.1 allows to obtain (see 
(4. 10), (4. 11)) the following estimate: 

e^° ^=** uo{x,t) dr. (5.2) 

-oo 

The above explicit corrispondence between uq and u^ clearly shows the interaction 
between diffusion and wave propagation. 

For instance, when fi — ^ sin(^), one has the pure wave 

uo{x,t) = sm{ — ) sin(— ci), (5-3) 

while Ue, owing to (5.2), is given by 

ueix,t) = e-fr^'' uoix,t). (5.4) 

This means that the dissipation caused by euxxt is significant by slow times e t only; 
the terms depending on the fast time t/e are fully negligible for all t>e. 
As consequence, in the time interval (e, -) where et < 1, the pure wave (5.3) is 
propagated nearly undisturbed; when t > -, damped oscillations predominate. 

These aspects are generally valid, whatever the initial data may be, and for 
appropriate /. In fact, if one transforms by (4.10) the time-components G'o,„(i) = 
sin(^ ct) of Go, one has [14]: 



(T-t)^ . ,TT 



^2„2 



„ sm{— cnr) dr — e 21^ ^ sin(— cnt). (5.5) 

Therefore, the Fourier series of the function H which approximates G^ is given by 

9 °° 1 
H= — y^ — e^2(^) s* sin(— cni) sm(—nx)sm(—n^). (5.6) 

rTT -*- — ^ r) I It 
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This formula, together with (2.6), (4.11), makes explicit the asymptotic behavior 
of the solution u^ of the problem Ve related to arbitrary data /o, /i and /. As (5.6) 
shows, this behavior is like (5.4), with damped waves which become vanishing as 
from the instants t > 1/e. 

It must be remarked that the function H can be given the form: 
H ^H~ -H+, with 



i/^ = — 2J- e"5( — ) s* sin(— ^)cos[— (x±ct)]. (5.7) 



CTT -"^-^ n I I 
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Then, the basic variables which are typical of the diffusion and wave propagation 
are identified by the following 

Property 5.1 - The function H~ (or H^ ) is solution of the diffusion-wave equation 

|.„..... K:|.„..-™.). ,.8, 

Moreover, if one puts: y = x ± ct, 9 = ^t, Eq. (5.8) turns into the heat equation 
Vyy = vg, with the time -variable 9 given just by the slow time and the space-variable 
y related to traveling (or retrograde) waves. 

Finally, we observe that H^ is given by Jacobi Theta function 9^{y,9); in fact, 
for y = X + ct or y = X — ct, one has: 

according to well-known formulae related to the heat equation. 
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